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Abstract 

Hydrogenation was carried out at different temperatures, 
hydrogen pressures and catalyst loadings and the observ- 
ations were used to tune parameters of a kinetic model. The 
reaction products were analyzed by Wijj's method to find 
iodine value, which quantitatively yielded the degree of 
unsaturation; hence conversion. Impure soy lecithin was 
hydrogenated using 10% Pd/C catalyst in pressure reactor at 
comparably low pressure such as 1.5 Mpa until iodine 
number values drops in the range of 10-30 gh/100g from 
initial 105gh/100g. The best fits were obtained for a reaction 
that occurs in three stages: adsorption of lecithin onto the 
catalyst surface (rate controlling step), hydrogenation of 
adsorbed lecithin and desorption of the hydrogenated 
product. The order of the rate determining step with respect 
to lecithin and catalyst is 0.81 and 0.53 respectively while the 
activation energy is 3.0 kcal/mol. The order with respect to 
hydrogen pressure is 0.6. The enthalpy of hydrogenation 
and desorption was determined to be -18.61 kcal/mol and 
20.35 kcal/mol respectively. The speed of agitation was 
sufficiently high enough to eliminate diffusional resistance. 
Hydrogenated lecithin provides products with increased 
stability with respect to oxidation, mostly applied in 
cosmetics and pharmaceutical preparations. 
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Introduction 

Lecithin, a natural emulsifying agent; naturally 
occurring mixture of the esters of stearic, palmitic and 
oleic acids, linked to the choline ester of phosphoric 
acid, plays an important role in cell metabolism and 
has great potential for use as a surfactant in cosmetics. 
Its use in cosmetics is limited due of its ready thermal 
and oxidative degradation: impure soy lecithin 

develops unpleasant odour upon its storage. The main 
cause of the instability is the presence of double bonds 
in the fatty acid esters present in the lecithin molecule 
which causes lecithin to develop rancidity on exposure 
to atmosphere. The double bond can be removed by 


one of several techniques including hydrogenation [Li 
& Xu (2009)], hydroxylation, acylation [Resh et al. (1999)], 
sulfonation [March (1985)], etc. These techniques 
modify the physical as well as chemical properties of 
this compound. Of these, selective hydrogenation of 
the double bonds yields a product that retains most of 
the desirable properties of the lecithin but with 
improved stability to oxidation and thermal degradation. 

The quality and physical properties of the final 
product during hydrogenation are greatly affected by 
the number of double bonds present in the oil, i.e. the 
iodine value and the cis-trans - isomers of fatty acid 
[Hunter (2005)]. 

Cole (1959) carried out hydrogenation of lecithin at 
low pressure (1.2-2.7 bar) in the presence of 10% 
palladium on charcoal catalyst in a mixture of 
chlorinated solvents such as methylene chloride and 
chlorobenzene and the temperature range from 0 to 
80°C. After filtration of the catalyst, chlorinated solvent 
and obtained product are with 60-80% yield. They also 
tried various other solvents such as benzene, ethyl 
acetate, methylene chloride, chlorobenzene and 
chloroform in order to increase the hydrogenation rate 
and thus reducing the operating pressure. However, 
the toxicity of chlorinated solvents and carcinogenic 
aromatics like benzene makes the process 
environmentally unfriendly and unsafe. Additionally, 
the stringent limits of the lowest possible content of 
the toxic solvents in the product put excessive 
economic burdens on the process. 

Mounts et al. (1978) carried out hydrogenation of 
Soybean oil performed in an autoclave with copper- 
chromite catalyst which included a high-pressure scan 
3.5, 7.0, and 21 MPa, at selected temperatures 110- 
170°C, agitator speed 1200-1400 rpm and at catalyst 
concentrations (0.05, 0.1, 0.2. and 0.4% copper). The 
samples were filtered to remove residual catalyst prior 
to analysis. Though the copper-chromite is most 
selective for hydrogenating the linolenic acid in 
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soyabean oil and other vegetable oils, drawback of 
commercially available copper-chromite catalysts is 
that they are much less active at low hydrogen 
pressures and spent catalyst presents a disposal 
problem. This is especially a problem since the 
hydrogenation of the lecithin requires higher 
temperature and pressure as compared to triglycerides. 

In view of the above discussion and the existing 
literature on hydrogenation of Impure Soy lecithin, 
there is need to study the hydrogenation process using 
environmentally friendly solvents and lower pressures 
and temperatures of operations. In this work a 
hydrogenation of Impure Soy lecithin over Pd/C 
catalyst with ethanol-toluene mixture as solvent. The 
advantage of this process is that the solvent is 
relatively benign and the catalyst can be recovered and 
reused. 

Materials and Methods 

Materials 

Impure Soy Lecithin was obtained from M/s Sonic 
Biochem, Indore, India. The rest of material such as 10% 
Pd/C (Particle size-6Micrometer), ethanol, toluene, 
acetone, Wijjs reagent, and all commercial grade 
chemicals were procured from standard reagent 
suppliers. All the reagents were laboratory grade. 

Methods 

In solvent mixture of ethanol/toluene (in the ratio of 
7:3), 5gm of Impure Soy Lecithin was dispersed along 
with 10% Pd/C catalyst in an autoclave and schematic 
diagram shown in FIG. 1. The catalyst loading of 0.05- 
2.5% (based on lecithin content w/w) and the 
temperature range 40-80°C were employed. The 
nitrogen gas pressure has been used for purging the 
reactor. During each run, the hydrogen partial 
pressure was held constant; the effect of hydrogen 
partial pressure was investigated for pressures 
between 0.3 MPa to 1.5 MPa. In all cases, the reaction 
was allowed to proceed for four hours. It was observed 
that the reaction became slower with an increase in 
toluene content in the solvent but we have not 
investigated this effect in detail. 

During the four hours of each run, 5ml aliquots of the 
reaction mixture were withdrawn every 30 minutes. 
Each sample was cooled to room temperature and it 
was diluted with toluene and filtered through Buckner 
funnel in order to remove the catalyst. The catalyst 
could be reused with some processing but the 
regeneration and activation of the catalyst is not in the 


scope of this work. The toluene layer obtained had 
acetone added to it and the liquid was concentrated by 
evaporation under vacuum to get a solid/gummy mass. 
It was later dried in a vacuum dryer to get the final 
product. 



FIG. 1 SCHEMATIC DIAGRAM OF THE EXPERIMENTAL SETUP 
A-HYDROGEN CYLINDER, B- NITROGEN CYLINDER, C- 
SAMPLING PORT WITH, D-STIRRER DRIVE SYSTEM, E- 
MAGNETIC DRIVE, F-PRESSURE AND TEMPERATURE 
CONTROLLER, G-TURBINE TYPE IMPELLER, H- DISCHARGING 
PORT, I- CHARGING PORT. 

Measure of Degree of Unsaturation (Iodine Value) 

Degree of unsaturation in the product was obtained 
from the iodine value according to the Wijj method 
[Jeeffery et al. (1991)]. During each hydrogenation run 
samples were withdrawn at 30 minute intervals. 0.2 
grams of the sample were mixed with 20 ml Wijj 
solution and 10 ml chloroform. It was then left in the 
dark for 30 minutes. Thereafter, 20 ml of 20% potassium 
iodide solution and 100 ml of distilled water was 
added. This was then titrated against 0.1M sodium 
thiosulfate (VI) solution. The experimental measurement 
error is +/- 0.05% for analysis of iodine value. 



FIG. 2 HPLC OF UNHYDROGENATED IMPURE SOY LECITHIN 
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FIG. 3HPLC OF HYDROGENATED IMPURE SOY LECITHIN 



FIG.4 INITIAL RATES OF REACTION (IN WIJJ VALUE PER 
SECOND) PLOTTED AGAINST IMPELLER SPEED (IN RPM). THE 
VARIATION WITH IMPELLER SPEED IS ONLY 3% BETWEEN 
1200 AND 1300 RPM. CONSEQUENTLY, THE REACTION IS 
CONSIDERED IMPELLER INDEPENDENT AT 1200 RPM. 
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FIG. 5 ORDER WITH RESPECT TO S. PLOT OF LOGARITHM OF 

INITIAL RATE LOG (- — ) VS LOGARITHM OF INITIAL 
V dt ) 0 

CATALYST CONCENTRATION LOG[5] 0 .THE SLOPE OF THIS 
PLOT IS M, THE ORDER WITH RESPECT TO [S], m = 0.41 + 

0. 13 (95% CONFIDENCE INTERVAL). 

High Performance Liquid Chromatography (HPLC) 

The analysis of Impure Soy Lecithin and hydrogenated 


Soy Lecithin was carried out using HPLC [Zhao et al. 
(2010)]. HPLC system (Agilent Technologies, India 
1200 Series) consisted of Chem Station software, model 
G1322A degasser, model G1311A quaternary pump, 
auto sampler, model G1316A column oven and 
variable wavelength UV-visible detector. The 
separation was performed on Agilent Eclipse plus C18 
column (4.6 X 100mm, 3.5pm). The mobile phase 
consisted of acetonitrile, methanol, and water in the 
ratio of 60:30:10 v/v. All the mobile phase solvents 
used were of HPLC grade filtered through 0.22 pm 
filter paper and sonicated for 20 min to remove any 
dissolved gases to avoid air entrapment in the HPLC 
system. Isocratic system of mobile phase was used for 
the elution of Impure Soy Lecithin at a flow rate of 
0.3ml min- 1 at constant column temperature of 30 °C 
and a wavelength of 210 nm was used for detection. 
All the samples were prepared by dissolving 
unhydrogenated and hydrogenated samples in 
methanol to get lmg/ml concentration and filtered 
through 0.22 pm Pall filter. 10 pi of sample was 
injected into the column. FIG. 2 shows number of 
impurities in impure soy lecithin with the PC, PE, and 
PI (see nomenclature) peaks. Unhydrogented PC 
shows the area of 1911.98 mAU in FIG. 4. From FIG.5, 
it can be said that the PC, PE, and PI shows retention 
time of 3.201, 4.455 and 5.858 respectively and others 
are impurities and to be removed by separation and 
purification process. Hydrogenated PC shows the area 
of 2627.8mAU in FIG. 3. 

Results and Discussion 

Agitation Speed 

The hydrogenation of Impure Soy Lecithin is a 
heterogeneous system involving two miscible solvents 
namely ethanol and toluene and the solid catalyst 
palladium on carbon. The diffusional process (mass 
transfer) of Impure Soy Lecithin and molecular 
hydrogen getting transported to active catalyst center 
depends upon the intensity of turbulence in the liquid 
phase [List et al. (1981)]. In order to study the effect of 
agitation, experiments at different impeller speeds 
were conducted at 80°Cand the initial rates of reaction, 
measured in terms of rate of reduction of Iodine (Wijj) 
value were plotted against the impeller speed as 
shown in FIG. 4. It was observed that the reduction in 
Iodine Value (Wijj Value) decreases rapidly with an 
increase in the speed of agitation until 1200 rpm, 
thereafter, the effect of agitation was minimal. Hence, 
all subsequent reactions were carried out at an agitator 
speed of 1200 rpm. 
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Reaction Mechanism 

The data was fitted using a kinetic model as described 
below. 

1) Adsorption of Lecithin (A) onto Catalyst (S) to 
Yield Adsorbed Lecithin (AS): 

A + 5 -> AS 

This reaction is assumed to be the rate controlling 
step. The order of reaction with respect to A is n 
and with respect to S is m. Since the mechanism 
within the catalyst particle is not studied in detail, 
the model is made with the "fudge" parameters m 
and n. In subsequent analysis, the step 1 will be 
further resolved to yield a sequence of elementary 
steps each with integral order of reactant. The rate 
of reaction with respect to A can be written as 
follows: 

77 = ~k A [A] n [S] m (1) 

Here, the kinetic constant is modeled using the 
Arrhenius equation as: 

k A = k Ao exp (- (2) 

Where k Ao and E A are constants to be fitted. 

2) Hydrogenation of Adsorbed Lecithin to 
Hydrogenated Product (CS) Using Hydrogen Gas (H2) 
under Pressure: 


H2 T AS CS 


This reaction is assumed to proceed rapidly and 
reaches its equilibrium condition. The equilibrium 

constant (. K AS ) can be written as: 

[cs] 


K t 


AS 


Ph 2 [as] 


(3) 


K as can be represented as a function of temperature 
by: 


A G° as = —RTlnK AS (4) 

Where A G AS is the free energy change associated 
with the hydrogenation reaction averaged over the 
temperature range of 40°Cto 80°C, which can be 
further expressed as 

A GSs = AH° s - T AS AS (5) 

Here A H AS and A S AS are the enthalpy and entropy 
change, respectively, accompanying the hydrogenation 
reaction averaged over the temperature range of 


40°C to 80°C. 


3) The Last Step in the Reaction Scheme Is 
Desorption of the Adsorbed Hydrogenated Lecithin 
(CS) to Give the Hydrogenated Lecithin Product (C) 
and Regenerate the Catalyst (S): 

CS c + s 

This step proceeds rapidly so that the equilibrium 
may be assumed. The equilibrium constant ( K cs ) 


can be written as : 


K, 


cs 


[C][S] 

[CS] 


(6) 


K cs can be represented as a function of temperature 
by: 


A G° s = -RTlnK cs (7) 

Where A is the free energy change associated 
with the desorption reaction averaged over the 
temperature range of 40°Cto80°C, which can be 
further expressed as follows: 

AG C ° S = A H° cs - T AS c S (8) 

Here A//° s and AS° 5 are the enthalpy and entropy 
change, respectively, accompanying the desorption 
reaction averaged over the temperature range. 
Additionally, assuming no hydrogenated lecithin is 
initially present; we have the following mass 
balance equations: 

[A \ 0 = [A] + [AS] + [CS] + [C] (9) 

Here, [A ] 0 is the concentration of lecithin taken 
initially. 


For the catalyst mass balance: 

[S ] 0 = [5] + [AS] + [CS] (10) 


Hence, from (3), (6), (9) and (10) we can write the 
following expressions for [C] and [S]: 



D] 0 — D] 

L J - k[S]+1 

(ll) 

And, 

rsi = [sl ° 

LJJ k[C]+1 

(12) 

Where, 

k — — — ( 1 H — j 

(13) 


KCS \ Ph 2 K AS J 


Eliminating [S] from equations (ll)and (12) gives. 


[C] 2 +a[C] -/? = 0 

(14) 

Where, 


a = ;+Mo-([4-W) 

(15) 

And, p = Wo_W 

K 

(16) 


Hence, using the quadratic formula and noting that 
[C] > o, 

[C] = - g+V f^ (17) 

The rate of reaction can be obtained by substituting 
equations (17) and (12) in equation (1). 


Fitting Data 

For each run, the concentration of lecithin in terms of 
Iodine Value was obtained at the interval of 30 
minutes. The results were used to tune the adjustable 
parameters in the reaction mechanism. To determine 
the order of reaction with respect to S, the logarithm of 
initial rates of reaction equation 1 was plotted against 
logarithm of concentration i.e. In vs ln[S] 0 for 
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temperature of 80°C and hydrogen pressure ( P Hz ) of 1.5 
MPa. The obtained plot gives a straight line, whose 
slope is m, the order of reaction with respect to [S], as 
shown in FIG.5. The value of m obtained from FIG.5 
was 0.41 + 0.13. 

Now, the order with respect to [A] is not straight- 
forward to determine since the initial concentrations of 
[A] are all the same for all our runs. The order with 
respect to hydrogen was similarly obtained as 0.6 (see 
figure 6). 

0.2 
o 

S 0 

T3 
< 

§ -0-2 
o 


TEMPERATURE). THE SLOPE OF THIS PLOT IS -E/JR, WHERE Ea 
IS THE ACTIVATION ENERGY FOR THE CONTROLLING STEP 
AND R IS THE UNIVERSAL GAS CONSTANT. HENCE— = 

R 

1388 ±205 K (95% confidence interval)I.E. £^ = 2.76 + 0.41 

Now, the rest of the parameters are not amenable to 
such exercise because the differential equation [1] is 
not tractable to solve analytically. Hence, we have 
solved it numerically using a Levenberg-Marquadt 
least-squares algorithm via the SciPy library of 
numerical tools [Jones et al. (2011)] in Python 
2. 6. However, we will compare the value of m and E A 
obtained from the numerical fit of all data with the 
value obtained from the analysis of FIG. 7 and as a test 
of the validity of the fitting. All graphs were plotted 
using the Matplotlib library of Python 2.6 [Hunter 
(2007)]. The result of the full numerical fit appears in 
Table 1. 


- 0.4 


log[P H2 ]o 

FIG. 6 ORDER WITH RESPECT TO HYDROGEN PRESSURE. PLOT 
OF LOGARITHM OF INITIAL RATE log (- VS LOGARITHM 
OF HYDROGEN PRESSURE log[P„ 2 ] THE SLOPE OF THIS PLOT 
IS THE ORDER WITH RESPECT TO HYDROGEN. 

There is another parameter that we can readily obtain 
from our preliminary data; and that is the activation 
energy E A/ from the Arrhenius plot (shown in FIG. 7) 
for initial rates of reaction for P Hz = 1.5 MPa and 
[S] o = 0.2 g-catalyst/lOOg-reaction mixture at 

temperatures between 40°C and 80°C.The slope of this 
plot is, —E a /R, where, R is the universal gas constant. 
Hence, the value obtained E A is 2.76 + 0.41^^ (95% 
confidence interval using student-t test). 
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FIG. 7 ACTIVATION ENERGY Ea. PLOT OF LOGARITHM OF 
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TABLEl VARIOUS PARAMETER FITTED VALUES OF KINETIC CONSTANTS. 


Fitted 

Parameter 

Optimized 

Value 

Unit 

Variation 

o 

< 

X 

0.1227 

mjj units) 1 - 71 ( 9 ~ cat )~ m 

K y VlOO g-mixj 




S 





2.76 ± 

Eao 

3.0 

kcal/mol 

0.41 (from 




FIG. 7) 

n 

0.81 






0.41 ± 

m 

0.53 


0.1 3 (from 




FIG. 5) 

AH» S 

-18.61 

kcal/mol 


AS£ s 

-0.0886 

kcal/(mol.K) 


AH» S 

20.35 

kcal/mol 


AS? s 

0.0985 

kcal/(mol.K) 



Discussion 

It was clear that the value of m and E A obtained from 
the full numerical fit is within the error of those 
obtained from linear regression, which validated our 
assumption for the adsorption of lecithin as the rate 
controlling step. The enthalpy of hydrogenation of 
lecithin appears to be -18.6 kcal/mol. Normally, 
enthalpies of hydrogenation of double bonds are 
around -30.0 kcal/mol [Hao et al. (2010), Timothy et al. 
(1988)]. The difference of about 10 kcal/mol can be 
straight forwardly attributed to the fact that the 
hydrogenation is not of the free lecithin but the 
adsorbed moiety and part of the bond-breaking may 
have occurred on the catalyst surface. The order of 
reaction with respect to lecithin was 0.81 with respect 
to the catalyst 0.53 and with respect to hydrogen 0.6. 
The reaction mechanism proposed was therefore not 
the whole story and we will subsequently work on 
elucidating a fuller mechanism by conducting 
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experiments with different catalyst morphologies. The 
activation energy of 2.76 kcal/mol was low compared 
to typical values greater than 8 kcal/mol. However, 
previous work [Hao et al. (2010)] has shown that 
activation energy for double bond hydrogenations 
may be lower than 1 kcal/mol. 

Conclusions 

Hydrogenation of Impure soy lecithin was carried out 
using solution of ethanol, toluene and palladium 
carbon as catalyst. The Pd/C catalyst was found to be a 
suitable catalyst for the hydrogenation of impure soy 
lecithin. The adsorption of lecithin onto the catalyst 
surface was determined to be the rate controlling step. 
The order of reaction with respect to lecithin 
concentration was 0.81 while that with respect to 
catalyst concentration was 0.53 and with respect to 
hydrogen was 0.6. The next step in the process 
elucidated the reaction steps that led to the fractional 
orders. Also, overall heat of reaction has been yet to be 
obtained. Once this parameter is obtained, the model 
may be used directly in reactor design. The kinetics 
gleaned from this work can be further used to design 
commercial scale reactors for mass production of this 
product using our catalyst. A more general model will 
be put forward in subsequent communications. 

Nomenclature 

Phosphatidylcholine (PC), Phosphatidyletanolamine 
(PE), Phosphatidylinositol (PI), Palladium carbon 
(Pd/C), High Performance Liquid Chromatography 
(HPLC), Activation Energy (EA),Milli absorbance 
(mAU), Rotation per minute (rpm). Hydrogen pressure 
(P H2 ), Mega Pascal(MPa). 
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